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Abstract 

Sun like stars are born from the collapse of fragment 
of molecular clouds. During the first embedded phase (the 
so called class 0/1), the protostar is surrounded by a col- 
lapsing envelope, whose physical and chemical structure 
would set up the initial conditions of the proto-stellar 
disk which may eventually form planets. In this contri- 
bution, we show how spectroscopic far infrared to mil- 
limeter observations can be a powerful diagnostic to de- 
rive the chemical and physical structure of the inner- 
most regions of envelopes collapsing around solar type 
protostars, where disk are believe to form. Our claim 
is based on the simultaneous modeling of the dynam- 
ics, chemistry and radiative transfer in a protostellar en- 
velope. Here we present a study of a class protostar 
NGC1333-IRAS4, obtained using ISO-LWS observations 
of water far-infrared lines. We show the derived physical 
and chemical structure of the protostar, namely the den- 
sity and temperature profiles, as well as the water abun- 
dance across the envelope. Besides, the mass of the cen- 
tral protostar and its accretion rate can be constrained by 
these observations and modeling, yielding an alternative 
method to measure these two poorly know parameters. 
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1. Introduction 

NGC1333-IRAS4 is a binary protostellar system, located 
in the south part of the Perseus cloud. The relatively 
large mass of the surrounding envelopes, together with 
the relative submillimeter to bolometric luminosity ratio 
suggest that both source are highly embedded and prob- 
ably very young. They have been classified as class pro- 
tostars (Andre et al. 1993). IRAS4 is associated with two 
molecular outflows, observed in the CO, CS (Blake et al. 



1995) and millimeters transitions (Lefloch et al. 1998). 
The dynamical age of these outflow is a few thousand 
of years. In the past years, several studies have focused 
on the molecular emission of this source. In particular, 
Blake et al. carried out a study of several molecules in 
NGC1333-IRAS4. Their major results were large deple- 
tion factors of all molecules in the envelopes, and the 
presence of a region with an increased abundance of CS, 
SiO and CH 3 OH that the authors attribute to grain man- 
tle desorption caused by grain-grain collisions induced by 
the powerful outflow originating from this source. More 
recently, Di Francesco et al. (2001) reported the detec- 
tion of inverse P-Cygni profiles of several lines toward 
NGC1333-IRAS4A, providing the least ambiguous detec- 
tion of infall motion towards a protostar ever. Here we 
present a study of the far infrared water line spectrum 
of NGC1333-IRAS4 observed with the Long Wavelength 
Spectrometer (Clegg et al. 1996, herein after LWS) on 
board ISO (Kessler et al. 1996). Using an accurate model 
of the thermal emission of the envelope, we show how the 
physical and chemical structure of the protostar can be 
derived (see also Maret et al. 2002). 



2. Observations 

A full range spectral survey of (43-196 /im) of NGC1333- 
IRAS4 was performed using LWS in the low resolution 
(R ~ 200) grating mode. The observations were obtained 
on three positions, the first one beeing centered between 
NGC1333-IRAS4A and NGC1333-IRAS4B, and the two 
other were taken along the outflow powered by NGC1333- 
IRAS4, NW and SW of it respectively. The Fig. 1 shows 
the spectra observed on the tree positions. We detected 
fourteen H2O lines, seven CO lines, three OH lines, the 
two [OI] and one [CII] fine structure line towards the cen- 
tral position. No significant H2O, CO and OH emission 
was detected on the NE and SW positions, whereas the 
[OI] and [CII] lines were detected towards the three po- 
sition at the same intensity 
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Figure 1. ISO-LWS spectra obtained towards NGC1333- 
IRAS4A on the central (top line), NE (middle line) and 
SW (bottom lines) positions. 



3. Line modeling 

The observed emission observed towards NGC1333-IRAS4 
can have a priori three different origin: the outflow pow- 
ered by this source, the PDR at the surface of the cloud, 
or the thermal emission of the protostellar envelopes. Ori- 
gin of molecular emission can be usually determined if a 
sufficient enough spatial or spectral resolution is achiev- 
able. In the case of ISO-LWS observations, all the gas 
components are mixed in the 80" beam of ISO. More- 
over, the relatively low spectral resolution does not al- 
low to spectrally disentangle these components. The out- 
flows powered by protostars has been invocated in the 
literature as a possible origin of the water line emission 
observed towards them (e.g. Giannini et al. 2001), but 
the absence of detection of water emission on the outflow 
lobes does not play in favor of this hypothesis. Moreover, 
Ceccarelli et al. (1999) noted the absence of the correla- 



tion between the H 2 and SiO emission 1 toward seven 
protostars, whereas they noted a correlation between the 
1.3 mm envelope emission. This suggest that the FIR wa- 
ter line emission can be due to the thermal emission of 
the envelopes surrounding the two protostars. 
In order to test this hypothesis, we modeled the wa- 
ter emission by an accurate model which simultaneously 
compute the thermal balance, chemistry and radiative 
transfer in the envelope (Ceccarelli, Hollenbach & Tie- 
lens 1996). The model has four free parameters: the cen- 
tral mass of the protostar M* , its accretion rate M , and 
the water abundance across the envelope. In order to take 
into account the increase of water abundance in the inner 
parts of the envelope - due to the evaporation of grain 
mantle when the gas reach the sublimation temperature 
of water ices - we used a different abundance in the in- 
ner (Xi n ) and outer parts (X out ) of the envelope. We 
run several models, varying the central mass between 0.3 
and 1 Mq , the accretion rate between 10~ 5 and 10~ 4 
Mq yr" 1 , and the water abundance in the outer and in- 
ner parts between 10~ 7 -10~ 6 and 10 _6 -10~ 5 respectively. 
The best fit parameters were obtained by minimizing the 
X 2 between the observations and the model results. The 
minimum \ 2 value is obtained for M* = 0.5 Mq , M = 5 
x 10~ 5 M yr" 1 , X out = 5 x 10~ 7 and X in = 5x 10~ 6 . 
The model reproduce relatively well the observed water 
emission, with the exception of two lines which seems to 
be underestimated by our model. 

4. Discussion 

The first result of this study is the capability of to repro- 
duce the observed water line emission with an envelope 
thermal emission model. Our model also allowed us to de- 
rive the density and temperature profile in the envelope. 
The Fig. 2 show the density and temperature profile de- 
rived. The observations are consistent with a Shu (1977) 
"inside-out" density profile, with a transition region start- 
ing at 1500 au. This profile is in good agreement with the 
one derived by J0rgensen et al. (2002), using continuum 
observations. The dust temperature closely track the dust 
one, with the exception of a small region at R <~ 50 au. 
This is causes at the radius where the temperature reach 
100 K, and the grain mantle evaporates, releasing large 
amount of water in the gas phase (about a factor ten 
more), increasing the cooling efficiency of the gas. This 

1 SiO is usually associated with the outflow strong shocks 
(e.g. Bachiller & Perez Gutierrez 1997) and it is believed to 
be a product of grain mantle desorption (Caselli et al. 1997; 
Schilke et al. 1997) 
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Figure 2. Density and temperature derived from the 
model. The solid line refer to the density profile in the 
envelope, the dashed line refers to the dust temperature, 
while the dotted line refers to the gas temperature. 



5. Conclusion 

Using a accurate model of the thermal emission of a 
protostellar envelope, we successfully modeled the water 
emission observed with ISO-LWS. We derived density and 
temperature profile in the envelopes, as well as the water 
abundance. The water abundance is enhanced by a factor 
ten in the inner parts, in an hot-core like region. We also 
derived the central mass of the protostar, which gives us 
an estimate of its age, if the accretion rate remains con- 
stant. Both mass, accretion rate, and water abundances 
compares well with recent studies. This work also em- 
phasis the need of ground based millimeter observation, 
where higher spatial and spectral resolutions are achiev- 
able. From this point of view, H 2 CO and CH 3 OH are 
of particular interest because of their large abundances 
in grain mantle. H 2 CO and CH 3 OH abundances should 
therefore be increased in the central parts of envelopes. 
This is indeed confirmed by very recent observations of 
these molecules towards several class protostars (Maret 
et al. in prep.). 



hot core like region, similar to the ones in massive star 
forming region - here indirectly probed by the jump in 
the water abundance needed to reproduce the observed 
emission - as been already observed on IRAS16293-2422, 
where a similar study was performed (Ceccarelli et al. 
2000a, 2000b). This claimed was confirmed by a recent 
study by Schoier et al. (2002). Furthermore, the central 
mass of the protostar is also constrained by our model, 
giving an alternative measure to this poorly known pa- 
rameters. We derived a central mass of 0.5 Mq and an 
accretion rate of 5 x 10~ 5 Mq yr _1 . This gives an age of 
10 4 year, close to the dynamical age of the outflow. This 
values compares well with the values derived recently by 
Di Francesco et al. The water abundance we found for the 
outer envelope is in good agreement with recent SWAS 
observations of the Perseus cloud (Bergin et al. 2002). As 
inner water abundances is dominated by the grain man- 
tle evaporation, the derived inner abundance should be 
a indirect measure of the water abundance in the grain 
mantle, the abundance we derive is somewhat lower than 
the water ice abundance observed towards massive star 
forming region (Tielens et al. 1984). Finally the abun- 
dances derived for NGC1333-IRAS4 arc in good agree- 
ment with the one derived by Ceccarelli et al (2000a) for 
IRAS16293-2422. This study show that NGC1333-IRAS4 
is probably younger than IRAS16293-2422. 
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